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A CLASS OF MULTIFRACTAL PROCESSES CONSTRUCTED 
USING AN EMBEDDED BRANCHING PROCESS 

By Geoffrey Decrouez and Owen Dafydd Jones 

University of Melbourne 

We present a new class of multifractal process on R, constructed 
using an embedded branching process. The construction makes use 
of known results on multitype branching random walks, and along 
the way constructs cascade measures on the boundaries of multitype 
Galton-Watson trees. Our class of processes includes Brownian mo- 
tion subjected to a continuous multifractal time-change. 

In addition, if we observe our process at a fixed spatial resolution, 
then we can obtain a finite Markov representation of it, which we can 
use for on-line simulation. That is, given only the Markov represen- 
tation at step n, we can generate step n + 1 in 0(log?i) operations. 
Detailed pseudo-code for this algorithm is provided. 

1. Introduction. Information about the local fluctuations of a process X 
can be obtained using the local exponent hx(t), defined as [37] 

hx(t) :=liminf log sup \X(u)—X(t)\. 

e^O loge | u _ t |< E 

When hx(t) is constant all along the sample path with probability 1, X is 
said to be monofractal. In contrast, we can consider a class of processes whose 
exponents behave erratically with time: each interval of positive length ex- 
hibits a range of different exponents. For such processes, it is, in practice, 
impossible to estimate hx(t) for all t, due to the finite precision of the data. 
Instead, we use the Hausdorff spectrum D(h), a global description of its 
local fluctuations. D(h) is defined as the Hausdorff dimension of the set of 
points with a given exponent h. For monofractal processes, D(h) degener- 
ates to a single point at some h = H [so D(H) = 1, and the convention is to 
set D{h) = — oo for H\. When the spectrum is nontrivial for a range of 
values of h, the process is said to be multifractal. 
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The term multifractal is also well defined for measures. Let B(x,r) be a 
ball centered at x € M. n with radius r. The local dimension of a finite measure 
jj, at x G M" is defined as 



The Hausdorff spectrum D(a) of a measure at scale a is then defined as 
the Hausdorff dimension of the set of points with a given local dimension a. 
Measures for which the Hausdorff spectrum does not degenerate to a point 
are called multifractal measures. Constructions of multifractal measures date 
back to the m-ary cascades of Mandelbrot [28], and the multifractal spec- 
trum of such measures can be found in, for example, [37]. 

A positive nondecreasing multifractal process can be obtained by inte- 
grating a multifractal measure. Other processes with nontrivial multifractal 
structure can be obtained by using the integrated measure as a multifractal 
time change, applied to monofractal processes such as fractional Brown- 
ian motion. This is the basis of models such as infinitely divisible cascades 



Multifractals have a wide range of applications. For example, the rich 
structure of network traffic exhibits multifractal patterns [1], as does the 
stock market [29, 30]. Other applications include turbulence [39], seismology 
[15, 40] and imaging [38], to cite but a few. 

On-line simulation of multifractal processes is in general difficult, because 
their correlations typically decay slowly, meaning that to simulate X(n + 1) 
one requires X(l), . . . , X(n). This is the same problem faced when simulating 
fractional Brownian motion, where to simulate X(n+ 1) one needs the whole 
covariance matrix of X(l), . . . ,X(n + 1). Some simple monofractal processes 
avoid this problem, for example, a-stable or M/G/oo processes [11], but 
it remains a real problem to find flexible multifractal models that can be 
quickly simulated. 

We propose a new class of multifractal processes, called Multifractal Em- 
bedded Branching Process (MEBP) processes, which can be efficiently sim- 
ulated on-line. MEBP are defined using the crossing tree, an ad-hoc space- 
time description of the process, and are such that the spatial component of 
their crossing tree is a Galton- Watson branching process. For any suitable 
branching process, there is a family of processes — identical up to a con- 
tinuous time change — for which the spatial component of the crossing tree 
coincides with the branching process. We identify one of these as the Canon- 
ical Embedded Branching Process (CEBP), and then construct MEBP from 
it using a multifractal time change. To allow on-line simulation of the pro- 
cess, the time change is constructed from a multiplicative cascade on the 
crossing tree. The simulation algorithm presents nice features since it only 
requires 0(n\ogn) operations and O(logn) storage to generate n steps, and 
can generate a new step on demand. 




[3, 5, 10]. 
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To construct the time change we use here, we start by constructing a 
multiplicative cascade on a multitype Galton-Watson tree. The cascade de- 
fines a measure on the boundary of the tree, whose existence follows from 
known results for multitype branching random walks. (See, e.g., [23] for 
the single-type case.) To map the cascade measure onto R+, we use the so 
called "branching measure" on the tree, in contrast to the way this is usu- 
ally done, using a "splitting measure." See Section 3 for details and further 
background. 

The MEBP processes constructed here include a couple of special cases 
of interest. We can represent Brownian motion as a CEBP, thus MEBP 
processes include a subclass of multifractal time changed Brownian motions. 
Such models are of particular interest in finance [29, 30]. In the special 
case when the number of subcrossings is constant and equal to two (for the 
definition see Section 2), the CEBP degenerates to a straight line, and the 
time change is just the well-known binary cascade (see, e.g., [4, 20, 32] and 
references therein). 

Although we do show that MEBP possess a form of discrete multifractal 
scaling [see the discussion following equation (3.5)], the multifractal nature 
of MEBP processes is not studied in this paper. We refer the reader to a 
coming paper for a full study of the multifactal spectrum of MEBP [13]. 
In particular, it can be shown that CEBP processes are monofractal, and 
that the multifractal formalism holds for MEBP processes, with a nontrivial 
spectrum. The monofractal nature of CEBP processes, together with an 
upper bound of the spectrum of MEBP, was derived in the Ph.D. thesis of 
the first author [12]. 

The paper is organized as follows. First we recall the definition of the 
crossing tree and then construct the CEBP process. We then construct 
MEBP processes and give conditions for continuity. Finally we provide an 
efficient on-line algorithm for simulating MEBP processes. An implementa- 
tion of the algorithm is available from the second author's website [16]. 

2. CEBP and the crossing tree. Let X : R + — >• R be a continuous process, 
with X(0) = 0. For n £ Z we define level n passage times TJ? by putting 
T n = and 

Jt +1 = inf{t > I%\X(t) € 2 n Z, X(t) + X(I?)}. 
The kth level n (equivalently scale 2 n ) crossing is the sample path from 
to T k ■ 

C%:={{t : X{t))\T%_ x <t<1t}. 

When passing from a coarse scale to a finer one, we decompose each level 
n crossing into a sequence of level n — 1 crossings. To define the crossing 
tree, we associate nodes with crossings, and the children of a node are its 



4 G. DECROUEZ AND O. D. JONES 




100 200 300 400 500 600 700 800 900 1000 




100 200 300 400 500 600 700 800 900 1000 

Fig. 1. A section of sample path and levels 3, 4 and 5 of its crossing tree. In the top frame 
we have joined the points T£ at each level, and in the bottom frame we have identified the 
kth level n crossing with the point (2 n ,T^ l _ 1 ) and linked each crossing to its subcrossings. 

subcrossings. The crossing tree is illustrated in Figure 1, where the level 3, 
4 and 5 crossings of a given sample path are shown. 

The crossing tree is an efficient way of representing a self-similar signal, 
and can also be used for inference. In [18] the crossing tree is used to test for 
self-similarity and to obtain an asymptotically consistent estimator of the 
Hurst index of a self-similar process with stationary increments, and in [19] 
it is used to test for stationarity. 

In addition to indexing crossings be their level and position within each 
level, we will also use a tree indexing scheme. Let be the root of the tree, 
representing the first level crossing. The first generation of children (which 
are level —1 crossings, of size 1/2) are labeled by i, 1 < i < Z , where Z& 
is the number of children of 0. The second generation (which are level —2 
crossings, of size 1/4) are then labeled ij, 1 < j < Zi, where is the number 
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of children of i. More generally, a node is an element of U = Un>o^ n an< ^ 
a branch is a couple (u,uj) where u£U and j G N. The length of a node 
i = i%, . . . ,i n is |i| = n, and the kth element is i[k] = ik- If |i| > n, i\ n is the 
curtailment of i after n terms. Conventionally \0\ = 0, and i|o = 0. A tree 
T is a set of nodes, that is, a subset of U, such that: 

• 0ST; 

• if a node i belongs to the tree, then every ancestor node k < |i|, belongs 
to the tree; 

• if u E T, then uj € T for j = 1, . . . , Z u and uj ^ T for j > Z u , where Z u 
is the number of children of u. 

Let T n be the nth generation of the tree, that is, the set of nodes of 
length n. (These are level — n crossings, of size 2~ n .) Define Ti = {j £ T||j| > 
|i| and j||;| =i}. The boundary of the tree is given by <9Y = {i £ N N |Vn > 
0,i| n € T}. Let be the position of node i within generation |i|, so that 

crossing i is just C J;]- The nodes to the left and right of i, namely C^^' ; 1 j_ 1 
and C^jj +1 , will be denoted i— and i+. In general, when we have quanti- 
ties associated with crossings, we will use tree indexing and level/position 
indexing interchangeably. So Z\ = Z^Yl , T\ = , etc. At present our tree 

indexing only applies to crossings contained within the first level crossing; 
however, in Section 3.3 we will extend this notation to the whole tree. 

Let € {+, — } be the orientation of + for up and — for down, and 
let A 1 ^ be the vector given by the orientations of the subcrossings of . Let 
= T£ — be the duration of C£. Clearly, to reconstruct the process 
we only need ct£ and for all n and k. The a£ encode the spatial behavior 
of the process, and the D"^ the temporal behavior. Our definition of an EBP 
is concerned with the spatial component only. 

Definition 2.1. A continuous process X with X(0) = is called an 
Embedded Branching Process (EBP) process if for any fixed n, conditioned 
on the crossing orientations the random variables A^, are all mutually 
independent, and A^ is conditionally independent of all A™ for m > n. In 
addition we require that {A 1 ^ = i} are identically distributed, for £ = +,—. 

That is, an EBP process is such that if we take any given crossing, then 
count the orientations of its subcrossings at successively finer scales, we get a 
(supercritical) two-type Galton-Watson process, where the types correspond 
to the orientations. 

Subcrossing orientations have a particular structure. A level n up cross- 
ing is from k2 n to (k + l)2 n , a down crossing is from k2 n to (k — l)2 n . 
The level n — 1 subcrossings that make up a level n parent crossing consist 
of excursions (up-down and down-up pairs) followed by a direct crossing 
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(down-down or up-up pairs), whose direction depends on the parent cross- 
ing: if the parent crossing is up, then the subcrossings end up-up; otherwise, 
they end down-down. Let be the length of that is, the number of 
subcrossings of C2- The number of up and down subcrossings will be written 
and Zu~, respectively. Clearly, each of the Z£ — 2 first entries of A% 
comes in pairs, each pair being up-down or down-up. The last two compo- 
nents are either the pair up-up or down-down, depending on at- Thus, given 
= +, we must have Z^ + = \Z"^ + 1 and Z^~ = \Z^ — 1, and conversely 
given a£ = — . 

Let A be the space of possible orientations. That is, a 6 A consists of some 

number of pairs, H — or — + , then a single pair +- 1- or . Given an EBP 

process, for the offspring type distributions we write p\{a) = P(^4J? = a\a^ = 
+) andp^(a) = P(AjJ = a\a% = -), for a E A. Let fi+ =E(Z"|a£ =+),/*- = 
E(Z^|a£ = — ) and \i = \{^ + + A* - ), then the mean offspring matrix is given 
by 

(z£+|< = +) (^rK = +)W^ + + i 
WK = -) (« = -)J V^-i + 

To proceed we need to make some assumptions about p A . 

Assumption 2.1. ^ + ,^ > 2 and E(Z™logZ™|a£ = j) < oo for i,j = ±. 

The first of these assumptions ensures that M is strictly positive with 
dominant eigenvalue \i > 2, and corresponding left eigenvector (ojo)- The 
corresponding right eigenvector is — 2)/( / u — 2),(/i~ — 2)/(// — 2)) T . 

The second assumption is the usual condition for the normed limit of a 
supercritical Galton-Watson process to be nontrivial. 

Theorem 2.1. For any offspring orientation distributions p A satisfying 
Assumption 2.1, there exists a corresponding continuous EBP process X 
defined on R + . 

Proof. A version of this result can be found as Theorem 1 in [17], for 
particular orientation distributions. 

Step 1. We initially construct a single crossing from to 1, with support 
[0,7?]. In step 2 we will extend the range to R and the support to [0, oo). 
X is obtained as the limit as n — > +oo, of a sequence of random walks X~ n 
with steps of size 2~ n and duration [T n . Put X°(0) = and X°(l) = 1, 
so that the coarsest scale is n = 0. Given X~ n we construct X~( n+1 ) by 
replacing the feth step of X~ n by a sequence of Z^ n steps of size 2~( n+1 ) and 
duration fi~^ n+1 \ If a^ n = i, then the orientations A^ n of the subcrossings 
are distributed according to p A . For a given n the A^ n are all mutually 
independent, and, given a^ n , A^ n is conditionally independent of all Aj m , 
for — m > —n. 
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Denote the (random) time that X n hits 1 by 

T 1 °'~ n = inf{t|X- n (t) = l}. 

We define X~ n {t) for all t £ [0, 7^ ,-n ] by linear interpolation, and set X~ n {t) = 
1 for all t > T^'~ n . The interpolated X~ n have continuous sample paths, and 
we will show that they converge uniformly on any finite interval, from which 
the continuity of the limit process follows. For any m < n, let r r~ m ~ n = o 
and 

T" + ™'~ n = mt{t > T~ m >- n \X- n {t) G 2- m Z,X- n (t) + X~ r ' '(T'" 1 '- 71 )} . 

UX~ n (T~ m ~ n ) = 1, then set T fc " + ™' _n = oo. By construction X- n (T~ m ~ n ) = 
X~ m (fi~ m k), for all k and m<n. The duration of the A:th level — m crossing 
of X~ n is D7 m >- n = T - m >- n - T~ m r n . 




Fig. 2. Construction of X° , X^ 1 and X~ 2 , and the associated crossing tree (see the 
proof of Theorem 2.1). The subtree rooted at crossing C^ 1 (the 4th crossing of size 1/2) 
corresponds to the tree inside the dashed box. In the notation of the proof of Theorem 2.1, 
for this subtree we have m = 1. If we go down one level in the subtree, corresponding to 
level n = 2 of the original tree, then Sti 4(1) = 3 and SZi 4 (1) = 1 count the number of 
up and down crossings at level 1 of the subtree. Similarly, 4 (2) = 7 and SZi 4 (2) = 3 
count the number of up and down crossings at level 2 of the subtree, and so on. This figure 
also illustrates other notation used in the proof of Theorem 2.1. For example, one has 



1,-2 



2 , since for X 2 , the the first crossing time of size 2 1 corresponds to the 



fourth crossing time of size 2 
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A realization of X°, X~ 1 and X~ 2 is given in Figure 2, with the associated 
crossing tree. 

We use a branching process result to establish that the crossing durations 
converge. When we defined the crossing tree (see Figure 1) we started with 
a sample path and then defined generations of crossings: taking the first 
crossing of size 1 as the root (level or generation 0), its subcrossings of size 
1/2 form the second generation (or level), its subcrossings of size 1/4 form 
the third generation, and so on. Each crossing can be up or down, so our 
tree has two types of nodes. Here we are reversing that process. That is, we 
are growing a tree using a two-type Galton-Watson process, and from the 
tree, constructing a sample path. The offspring distributions for our tree are 
just p^. Given the tree at generation n, we get an approximate sample path 
by taking a sequence of up and down steps of size 2~ n and duration /x~ n , 
with directions taken from the node types of the tree. We need to show that 
the sequence of sample paths, obtained as n — > oo, converges. 

Consider the subtree descending from crossing C k m . Let S_ mk (n — m) 

and SZ mk (n — m) be the number of up and down crossings of size 2~ n 
which are descended from the kth crossing of size 2 _m ; then {(S^ mk (n — 
m),SZ m k( n ~ m ))}'n=m 1S a two- type Galton-Watson process. From Athreya 
and Ney [2], Section V.6, Theorems 1 and 2, we have that as n — > oo, 
fj m ~ n [SZ m k (n — m),SZ m k(. n ~ m )) converges almost surely and in mean to 
(i, l)W k ~ m , where W k m is strictly positive, continuous and ¥,{W k m \a k m = 
±) = (jj± - 2)/Gu - 2). Mor eover, the distribution of W k m depends only 
on a k m , and for any fixed m the W k m are all independent. Finally, since 
S - m ,k( n ~m) + SZ m , k (n -m) = fi n D~ m ~ n , we have 

D k m, ~ n -> V~ m W k m a.s. as n ^ oo. 

Accordingly, let T k m = £* =1 fi~ m Wf m = ]im n -^ 00 T~ m, ~ n . 

Take any e > 0, 5 > and T > 0. To establish the a.s. convergence of the 
processes X~ n , uniformly on compact intervals, we show that we can find a 
u so that with probability 1 — e, 

(2.1) \X~ r (t) - X~ s (t)\ <5 for all r,s >« and te[0,T]. 
Given t € [0,T], let k = k(n,t) be such that 

rp — n ^ ± rp — rl 

1 k-l- l< - 1 k ■ 
For any r,s >n, the triangle inequality yields 

\X~ r {t) - X~ s {t)\ < \X~ r {t) - X- r {T7 n '~ r )\ 

(2.2) 
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since X~ r {T~ n ~ r ) = X~ s {T~ n ~ s ) = X~ n {k^- n ). 

For any u > n let j = j(n, u) be the smallest j such that r p~ n ~ u > T. As 
u — > +00, j(n,u) — > j(n) < 00 a.s., so for any n we can choose Eq such that 

p( min fi~ n Wr n >e )>l-e, 

\i<j(n) / 

and u such that for all q>u, 

w(m^\T- n '- q -T- n \<e )>l-e, 

\i<](n) / 

which yields 

Pf max \T~ n '- q -Tr n \< min fi~ n Wr n ) >l-e. 

\i<j(n) i<j(n) / 

Thus, given n we can find u such that for all q>u, with probability at least 

1-8, 

T k-f q < 1 < T k+f q fOT a11 1 € [°' T ] • 

Now, since X^^f 11 ) = X~ n {{k - 2)n~ n ), X^(T^f q ) = X~ n ((k + 
l)[i~ n ), and in three steps X~ n can move at most distance 3 • 2~ n , we have 

\X-i(t) - X-i(T- n >- q )\ < 3 • 2" n . 

Choosing n large enough that 6 • 2~ ra < 5, we see that (2.1) follows from 

(2.2) . Sending 5 and e to shows that X~ n converges to some continuous 
limit process X uniformly on all closed intervals [0, T], with probability 1. 
By construction, the duration of crossing C^ n is pT n W^ n . 

Step 2. Clearly the construction above can be used to generate any cross- 
ing from to ±2 n . Thus, to extend our construction from a single crossing 
to a process X(t) defined for all t € R + , we proceed by constructing a nested 
sequence of processes {X^J-^Lq, such that X^ is a crossing from to ±2 n , 
and the first level n crossing of X^ n+1 ^ is precisely X^ n \ To make this work, 
we just need to specify P(X^ n ^(T") = 2 n ) in a consistent manner. 

Consider the orientation of the first crossing from to ±2 n for an EBP 
process. Let u = P(a? = +K +1 = +) and v = P(a? = +K +1 = -); then u 
and v are determined by , and 

(2.3) a n := P(q™ = +) = ua n+ i + v(l - a n+ i) = v + (u- v)a n+i . 

For (n, v) € [0, l] 2 \ {(1,0)}, we see that equation (2.3) has fixed point a = 
v/{l — u + v) € [0,1]. Moreover, the only doubly infinite sequence {a n }'^ = _ 00 
which satisfies (2.3) and remains in [0, 1] is given by a n = a for all n. Given 
this, it follows that a n = a, and thus from Bayes's theorem that P(a n+1 = 
+ | a J = +) = u and F(a n+1 = +K = -) = v. If (u,v) = (1,0), then any 
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a € [0, 1] is possible, but everything else goes through as before. In this case 
the a™ are all the same, but may be of either type. 

Construct as a crossing from to 1 with probability a [the fixed point 
of (2.3)], otherwise as a crossing from to —1. Then, given X^ n \ construct 
X(n+l) ag follows: first, put = + with probability u if a\ = +, with 

probability v otherwise; second, generate A™ +1 conditional on a™ +1 and a"; 
third, use as the first level n crossing of X^ n+1 ^; finally construct the 
remaining level n crossings conditional on cxgjCKg, . . . ,d r t n+1 . Write X for 

the limit of the X^ n \ To complete our construction we just need to check 
that the process X does not escape to ±oo in finite time. By construction, 
we have T? = ml{t\X(t) = ±2 n } = fj, n W?, where W? is strictly positive, 
continuous, and has a distribution depending only on the orientation a". 
Thus for any T > 0, P(Tf < T) -> as n -> oo. □ 

Theorem 2.2. Let X be the EBP constructed in Theorem 2.1; then, 
for each n, conditioned on the crossing orientations a 7 ?, the crossing du- 
rations are all mutually independent, and is conditionally inde- 
pendent of all AJ 1 for m>n. Also, E(D£|a£ = ±) = ^"(/i* - 2)/(/i - 2), 
and the distribution of fi~ n D r k l depends only on a^. Moreover, up to finite- 
dimensional distributions, X is the unique such EBP with offspring ori- 
entation distributions p^. That is, for any other EBP process Y with off- 
spring orientation distributions p^ and crossing durations as above, we have 

(X(h), . . .,X(t k )) = (Y(h), Y(t k )) forany0<t l <t 2 <---<t k . 

Accordingly, we call X the Canonical EBP (CEBP) process with these 
offspring distributions. 

We also observe that X is discrete scale-invariant: let H = log 2/ log /x; 
then for all c 6 {/i n ,n £ Z}, 

(2.4) X(t) f = c~ H X{ct), 

where = denotes equality for finite- dimensional distributions. H = log /x/ log 2 
is known as the Hurst index. 

Proof. We retain the notation of Theorem 2.1. 

For the process X, the dependence structure of the crossing durations is 
clear from the construction. 

To show uniqueness, let Y be some other EBP process with offspring 
orientation distributions p^ , and crossing durations satisfying the conditions 
of the theorem statement. We will make use of the same notation for the 
crossing times, durations, orientations, etc. of Y as for X, and rely on the 
context to distinguish them. 

For an EBP, the finite joint distributions of the orientations A 1 ^ are de- 
termined completely by p^, and thus are identical for X and Y. For the 
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crossing durations of Y, note that for any m<n and k, we have 

f (— 77i,— n,fc+l) 

(2.5) ^ m D- m = fx m - n V n D- n , 

j=C(-m,-n,k)+l 

where £(— m, — n, fc) is such that £(—m,—n,k) + 1 is the index of the first 
level — n subcrossing of C^ m . Thus by the strong law of large numbers, 
sending n oo, 

C(-m,~n,k+X) 



j'=f (— 77i,— n,fe)+l 



+ /z m -"5Z m , fc (n-m) 



n C(- m :- n : fc + l) 

S -m,k( n ~ m ) j=C(-m ,-n,k)+l 



A ^ fc - m M n E(Z)- n |a7 n = +) + ^ fe - m ^E( J D7™| a 7™ = _) 

= wr", 

where the distribution of Wj~ m is completely determined by pj^ , and thus is 
the same for X and Y. 

Once we have the crossing orientations and the assumed dependence struc- 
ture of the crossing durations, the crossing distributions (for up and down 
types) determine the joint distributions of the crossing times {TJ}}. Thus, 
for any n and k, {X(Tp)} k =0 and {Y(T™)} k =0 are identically distributed. 
Since any t can be bracketed by a sequence of hitting times, X and Y are 
identical up to finite-dimensional distributions. 

That X is discrete scale-invariant is a direct consequence of its construc- 
tion, since simultaneously scaling the state space by 2 k and time space by 
/j, k does not change the distribution of X. □ 

Remark 2.1. From [14] it is clear that Brownian motion is an example 
of a CEBP process, where the offspring of any crossing consist of a geometric 
(1/2) number of excursions, each up-down or down-up with equal proba- 
bility, followed by either an up-up or down-down direct crossing. That is, 

p!(---++)=Pa( ) = 2"(* +1 \ 

where • ■ • represents a combination of z pairs, each either -| — or — h It 
follows that P(Z£ = 2x) = 2~ x , independently of ajj. 

3. From CEBP to MEBP. In this section we construct Multifractal Em- 
bedded Branching processes (MEBP processes) as time changed CEBP pro- 
cesses. 
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Consider initially a single crossing of a CEBP X, from to ±1. We 
constructed X as the limit of a sequence of processes X~ n , which take 
steps of size 2~ n and duration p~ n . The crossing tree gives the number of 
subcrossings of each crossing. If we add a weight of l//i to each branch of 
the tree, then truncating the tree at level — n, the product of the weights 
down any line of descent is /x _n , which is the duration of any single crossing 
by X~ n . We generalize this by allowing the weights to be random, then 
defining the duration of a crossing to be the product of the random weights 
down the line of descent of the crossing. The resulting process, Y~ n say, can 
be viewed as a time-change of X~ n , where the time-change is obtained from 
a multiplicative cascade defined on a (two-type) Galton-Watson tree. 

As for CEBP, we will initially construct a single level crossing of an 
MEBP, then extend the construction to R+. We will retain the notation 
of Section 2, but note that we will prefer the tree indexing scheme to the 
level/position indexing scheme in what follows. In particular, the number 
of level —n up and down subcrossings of node i in level —m are denoted 
S^~(n — m) and S7~(n — m), and, under Assumption 2.1, the almost sure limit 
and mean limit of /i m_n (5 ; + (?i — m),S^(n — m)) is (^, The duration 

of crossing i of the CEBP process X is then ^i~ m W\. 

We assign weight R\{j) to the branch (i, ij). R\ := {R\{1), . . . , R\(Z{)) may 
depend on A\, but conditioned on a\ must be independent of other nodes 
j that are not descendants of i. For r 6 R^ , write F^ a (r) = ^(^(l) < 
r(l), . . . ,R (z) < r(z)\a = ±,A = a, \a\ = z) for the joint distribution of 
R , conditioned on the crossing orientations a. The weight attributed to 
node i is 

|i|-i 

Pi =n + !])• 

k=0 

That is, pi is the product of all weights on the line of descent from the root 
down to node i. We use the weights to define a measure, u, on the boundary 
of the crossing tree. The measure v on <9T is then mapped to a measure £ 
on R, with which we define a chronometer A4 (a nondecreasing process) by 
M(t) = C([0,i]). The MEBP process is then given by Y = XoM" 1 , where X 
is the CEBP. The crossing trees of X and Y have the same spatial structure, 
but have different crossing durations. In Figure 4 we plot a realization of an 
MEBP process and its associated CEBP. 

The literature on multiplicative cascades is rather extensive. For the ex- 
istence of limit random measures and the study of the properties of cer- 
tain martingales defined on m-ary trees, one can refer, for instance, to the 
works of Kahane and Peyriere [20], Barral [4], Liu and Rouault [25] and 
Peyriere [36]. For results on random cascades defined on Galton-Watson 
trees, see, for example, Liu [22, 23], and Burd and Waymire [9]. 
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To obtain the time-change process explicitly, the random measure defined 
on the boundary of the tree is mapped to M + then integrated. Note that 
this mapping, given explicitly in Section 3.2, differs from random partitions 
previously considered in the literature. The usual approach is to use a "split- 
ting measure" to map the boundary of the tree to [0,1], then use the density 
of the cascade measure with respect to the splitting measure; see, for exam- 
ple, [34, 35, 37]. Our approach can be thought of as using the "branching 
measure" instead of a splitting measure. A splitting measure is constructed 
by splitting the mass associated with a given node between its offspring, 
with no mass lost or gained. The branching measure allocates mass accord- 
ing to the number of offspring, and is only conserved in mean. We have 
taken the terminology of splitting and branching measures from [23] , Exam- 
ple 1.3. A multifractal study of the measure we construct on M + is given in 
a forth-coming paper [13]. 

3.1. The measure v. To construct u, we use a well-known correspon- 
dence between branching random walks and random cascades, in which the 
offspring of individual i have types given by A{ and displacements (relative to 
i) given by — log-Rj. For background on multitype branching random walks, 
we refer the reader to Kyprianou and Sani [21] and Biggins and Sani [8]. 

Suppose |i| = m and n>m. Define 



Wt(n-m,6)= V (pj/pO 



jGT n nTi,aj=± 



and for i,j 



m i , J (6) = E(W (l,9)\a = i) 



E 



( ]T R {k) e \a = i\ 



Let M{6) = (mij(9))ij = ±, and write mfj(9) for the entry of the nth 

power M n {9). Then it is straight forward to check that E(Wj (n — m, 0)\a\ = 

i) = m™~ m (6). If we take constant weights equal to 1/u, then W^(n — m, 1) = 

v 

fi m n S i , in the notation of Theorem 2.1. 

Let (i(6) be the largest eigenvalue of M(0). We make the following as- 
sumptions about R . 

Assumption 3.1. We suppose that < R < oo a.s., M(8) < oo in an 
open neighborhood of 1, /i(l) = 1 and u'(l) < 0. 

In the case where the distribution of R (and thus Z ) does not depend 
on a , we assume in addition that 

^i? (j)bg^ (i) <°o. 

j=l j=l J 
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In the case where there is dependence on the crossing orientation (type), 
we suppose that for some 5 > 1, fi(5) < 1 and 

\ 6 \ 




a = i I < oo for i = ±. 



Note that if the weights are finite and strictly positive, then M and // 
from the previous section are just M(0) and //(0), and from Assumption 2.1 
we get < M(9) for all > 0. In the case where R does not depend on 
a , the BRW simplifies to a single-type process, and the condition fj,(l) = 1 
simplifies to IEQJ^j Rmij)) = 1, which we recognize as a conservation of 
mass condition. 

Left and right eigenvectors corresponding to /i(l) will be denoted u = 
(u + ,u~) and v = (v + ,v~) T , normed so that u(l, 1) T = 1 and uv = 1. The 
following lemma is a direct consequence of Biggins and Kyprianou [7] , The- 
orem 7.1, and Biggins and Sani [8], Theorem 4. 

Lemma 3.1. Under Assumptions 2.1 and 3.1, (W; + (n — m, l),W;~(n — 
m, 1)) converges almost surely to uW;, for some random variable Wi suc/i 
that the distribution q/Wi depends only on a\, and E(Wi|ai = i) =v l . More- 
over, for each n, conditioned on the crossing orientations ai, i S T n , the Wi 
are mutually independent, and Wi is conditionally independent of (Aj,Rj) 
for |j| < |ij. For all nodes i, 

(3.1) W i = ^E i (i)W ii . 

j'=i 

Note that in the case where R does not depend on a , the right eigen- 
vector v = (1, 1) T . 

We can now define the measure v on <9Y. Recall Yj = {j € Yj|j| > |i| 
and j||i| = i}, so d1\ contains all the nodes on the boundary of the tree 
which have i as an ancestor. We define i^((9Yi) = p(W\. By Caratheodory's 
extension theorem, we can uniquely extend v to the sigma algebra generated 
by these cylinder sets. 

3.2. The measure £ and time change M. The measure £ is a mapping 
of v from <9Y to [0, W ] C M. By analogy with m-ary cascades, we call £ a 
Galton- Watson cascade measure on [0, W ]. 

As above, let T^ n denote the /cth level — n passage time of the CEBP 
process X, and put 
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Putting C({0}) = 0, this gives us (([0,TjT n ]) for all n,k > 0. For arbitrary 
t G (O,W ], let i G <9T be such that t G (T"" _ 1 ,T~" ] for all n > 0. 

Noting that T^j » = Tj| n is a nonincreasing sequence, we define £([0,t]) = 

iim n _, oo c([0,r- ( ^ ) ]). 

We can now define .M(t) = £([0,i]), and define the MEBP process Y (on 
[O,W ]) as 

y = ioM 4 . 

Here we take .A/f _1 (i) = inf{s : M(s) > i}, so that it is well defined, even if 
M. has jumps or flat spots. 

Put T k 7 n = M(T~ n ) = £* =1 />7 n Wr™. Then y(7^ n ) = X(T~ n ), so T k 7 n 
is the feth level —n crossing time for Y, and T> k n = p k n W k n the kth. level 
— n crossing duration. Note that if we take constant weights equal to 1/fJ,, 
then T~ n = T~ n and Y = X. 

Lemma 3.2. Under Assumptions 2.1 and 3.1, A4 and M^ 1 are contin- 
uous. That is, M has neither jumps nor flat spots. 

Proof. To show that A4 has no flat spots, it is enough to show that: 
(a) 

— -|) TT? Tl IP fx 

max /i W k — > asra— >-oo, 

(b) 

W k n > a.s. for each n, k > 0. 

Property (a) follows directly from Theorem 1 in [33], noting that under 
Assumption 2.1 E(Wj|ai = ±) < oo, so that xdF ± (x) is slowly varying, 
where i ?± (x) = ¥{W\ < x|aj = ±). This is equivalent to saying that the mea- 
sure P, defined on <9Y by u(dT{) = fi'^Wi, has no atoms. 

To show (b), let q± = P(W = O|a = ±), then note that since the weights 
R > 0, we have, from (3.1), that 

(3-2) Qi = fi(q+,Q-), 

where fa is the joint probability generating function of Z@ given a = 
i. [Note that q± = F(W = O|a = ±) satisfy the same equations.] Since 
(Z l \a = i) > 2 and F(Z± = 2\a = i) < 1, we have for G [0, l] 2 \ 

{(0,0), (1,1)}, < Thus the only solutions to (3.2) are (0,0) and 

(1, 1), and as E(W |a = ±) > 0, we get q± = 0. 

A4 is continuous (has no jumps) if £ has no atoms. That is, 

(a*) 

max p k n y\> k n as n -> oo, 
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(b*) 



W k n > a.s. for each n, k > 0. 



We prove (b*) in exactly the same way as (b). 

Property (a*) is equivalent to saying that v has no atoms. In the case 
where the distribution of R does not depend on a&, the BRW embedded 
in the crossing tree is effectively single-type, and (a*) is given by Liu and 
Rouault [24], Theorem 6. In the case where the distribution of R does 
depend on a , the approach of [24] generalizes only as far as the end of 
their Lemma 13, at which point we require, for some A < 1, 

(3.3) Ei/({i:pi| n > A n }) asn^oo. 

However, this can be shown using some recent results of Biggins [6], as we 
now demonstrate. 

In the notation of [6], consider a BRW with offspring types {cij} = {A (i)} 
and displacements {zi} = {log(R (i) / y)}, for some 7 > 1. Put 

\<k<Z ,A (k)=j ' 

(this is 771 j j in the notation of [6]). Then the matrix M{8) = (rhij(9))ij = ± 
has maximum "Perron-Frobenius" eigenvalue k(9) = p{9)^ 9 . From assump- 
tions 2.1 and 3.1 it is clear that for some 9 > 0, M(9) is finite, irreducible 
and primitive. 

(n) 

Let B\ be the rightmost particle of type % in generation n, that is, 
B±^ = max log p\ + n log 7. 

ieT n ,a ; =± 

Then Proposition 5.6 of [6] shows that 



B {n) 



n 

where K*(a) = snp d>0 {9a — k(8)} and T(k*) = sup{a : K*(a) < 0}. 

We have k(0) = /^(0), k(1) = 7, ft'(l) = /i'(l)7 + 1, and for 7 large enough, 
k(9) — >• 00 as 9 —> 00, faster than linear. T(k*) corresponds to the slope of the 
line that passes through the origin and is tangent to k* , from which it follows 
that r(«*) < 7 provided that k'(1) / 7, that is, provided / (7 — 
But //(l) < and 7 > 1 by assumption, so ^ 7, and we get 

log Pi a.s , ^ 
max > L(k ) — 7 < U. 

i6T n , ai =± n 

Equation (3.3) follows immediately, completing the proof of our lemma. □ 

3.3. Extending the construction to M+. We can extend Y from [0, W0] 
to in much the same way we extended the CEBP X, by constructing a 
sequence of nested processes Y^ n \ where y( n ) consists of a a single level n 
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crossing from to ±2 n , and the first level n crossing of y( n+1 ) is precisely 
YH As for the CEBP we need to specify P(y( n )(7^ n ) = 2 n ) in a consistent 
manner, but we also need to scale the first crossing. 

Construct Y (0) as a crossing from to 1 with probability a [the fixed 
point of (2.3)], otherwise as a crossing from to —1. Then, given 

Y (n) 

construct y( n+1 ) as follows: first, put a™ +1 = + with probability u if a™ = + 

and probability v otherwise; second, generate (A™ +1 ,R™ +1 ) conditional on 

< +1 and a?; third, scale the weights R™ +1 by fourth, use 

as the first level n crossing of Y ( n+1 ) ; finally, construct the remaining level n 

crossings conditional on , o;g, . . . , a™ +1 . Write Y for the limit of the Y^ n \ 

i 

When constructing y( n+1 ) we take independent processes, each con- 
structed like 7^, then scale the first by 1 = E% +1 (1) / (1) , the second 
by R™ +1 (2) / R™ +1 (1) , and so on, before stitching them together. When con- 
structing the second and subsequent level n crossings 

of y(n+l) 

we proceed 

exactly as for the construction of Y ^ , except for a spatial scaling of 2 n and 
a temporal scaling of ni=i V-^i(-Q) n °ting that the R\{^) are taken from 
the first level n crossing, and are thus independent of the second and subse- 
quent level n crossings. Thus with this construction, the process Y^ n \t) is 
distributed as 2 n Y^ (tpj~ n ), where p^ n is the weight given to the first level 
— n crossing of Y^ (a product of n weights, from level —1 to — n). 

To complete our construction, we just need to check that the process Y 
does not escape to ±oo in finite time. To see this note that the second level 
n crossing of y( n+1 ) is distributed as 

J R™ +1 (2) 



where, conditioned on its orientation, is equal in distribution to the 
level crossing of Y^°\ and is independent of i?J (1) for k = 1, . . . ,n + 1 
and of i?™ +1 (2). We have already seen that > almost surely, and by 
assumption, i?™ +1 (2) > 0, so it suffices to show that Ilfeii ^i (1) ~~ * almost 
surely as n — > oo. Given the orientations af, k = 1, . . . ,n + I, the weights 
Ri(l) are independent. The sequence of orientations {ceil^l! form a two- 
state (+ and — ) Markov chain, with transition matrix 

u 1 — u 
v 1 — v 

Thus the product R\(\)R\{\) ■ ■ ■ can be written as a product of independent 
random variables of the form 

U V 
C=J{A k J{B k , 

k=l k=l 
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where U ~ geom(u), V ~ geom(l - v), A k ~ (R\ {l)\a\ = +), B k ~ (i?i(l)| 
a\{l) = — ), and they are all independent. The product rifeii-^i(l) con ~ 
verges to zero if the sum X/fcii \ogR\ (1) diverges to — oo, which follows 
almost surely from the strong law of large numbers, provided ElogC = 
Y^Elog^4i + -ElogBi < (assuming u^l and v ^ 0). That is, the pro- 
cess Y is defined on IR + provided the following assumption holds. 

Assumption 3.2. If u = P(a? = +K +1 = +) + 1 and v = P(a? = 
+|a" +1 = — ) 7^ 0, then we suppose that 

_LE(logi??(l)K = +) + -E(logi^(l)|a? = -) < 0. 
1 — u v 

If u = 1, then we require E(logi?™(l)|a™ = +) < 0, and if v = 0, then we 
require E(logi^(l)K = -) < 0. 

To describe the crossing tree of the extended process Y, it is convenient to 
extend the tree-indexing notation introduced earlier. We do this by indexing 
nodes relative to a spine, defined by the first crossing at each level. For any 
node in the tree, we can trace its ancestry back to the spine. For any n let 
n : be the node on level n of the spine and T n . the tree descending from 
that node. Nodes in the tree T n . will be labeled n : i, where i is the node 
index relative to n : 0. Thus n : i is in level n — |i| of the crossing tree, and 
a crossing previously labeled i is now labeled : i. Note that this labeling is 
not unique, as n : i = (n + 1) : li. 

Write p n: i for the weight assigned to node n:i, which is given by 

|i|-l ,("A|i|)-l 

H R n . Alk (i[k + l]) / 11 i? (n _ fc):0 (l), n>0, 

k=0 ' k=0 

n^.i^fc + l]), n = 0, 

k=0 

|i|-i H-i 

Y[R n:i{k (i[k + l}) J] i^ fc!0 (l), n<0. 

k=0 k=0 

Here we have used the convention that nfc=o Xfc = 1, to deal with the case 
|i| =0. 

Let W n: i be branching random walk limit associated with crossing n:i; 
see Lemma 3.1. Then the duration of crossing n : i is 

(3.4) V n . A = p n:i W n:i . 

We summarize conditions for existence and continuity of Y in the theorem 
below. 
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Theorem 3.1. Suppose we are given subcrossing orientation distribu- 
tions p^ and weight distributions F R ,, satisfying Assumptions 2.1, 3.1 
and 3.2. Then there exists a continuous EBP process Y with subcrossing 

orientation distributions p^ and crossing durations T> n . -J= p n - iW n: i. 

For each n, conditioned on the crossing orientations a^, the random vari- 
ables Wj! are mutually independent, and W% is conditionally independent 
of all (A™,R™) for m>n. Also, E(W£|a£ = i) = v\ and the distribution of 
W% depends only on oUt. 

We call Y the multifractal embedded branching process (MEBP) defined 

by Pa andF %a 

As a corollary of our construction we also obtain a novel Galton- Watson 
cascade measure £ on R + , constructed by mapping the cascade measure u 
from the boundary of the (doubly infinite) tree to R+, using the measure u 
as a reference. [Where v is defined on <9Y by P(<9Y n: i) = 

Mandelbrot, Fisher and Calvet [31] described a class of multifractal pro- 
cesses such that 

Y(at) f = M{a)Y{t) and M(ab) = M 1 (a)M 2 {b), 

where M\ and M 2 are independent copies of M. Write A for M , and then 
we can re-express the scaling rule for Y as 

(3.5) Y(A{a)t) f =aY(t) and A{ab) = A l (a)A 2 (b), 

where A\ and A 2 are independent copies of A. When constructing our 
MEBP Y, we noted that Y^ n \t) is distributed as 2 n Y^(tp~ n ). More gen- 
erally we have y( m+n )(t) = 2 n Y^ m \tp^[ n ), so sending m — > oo we get, for 
n = 0,l,..., 

Y (t) i =2 n Y(tpi n ). 

This is close to the form (3.5) with A(2~ n ) = p^ n = ^lC 1 )- The de- 

ferences are that A(a) is only defined for a = 2 _n , n 6 Z+, and the product 

form A(ab) = Ai{a)A 2 {b) does not quite hold because of the dependence of 
R\ on the orientation a\. [In fact, the sequence {(— log p^ n , a^ n )} is Markov 
additive.] Nonetheless, we recognize that MEBP processes possess a form of 
discrete multifractal scaling. The full multifractal spectrum is obtained in a 
forthcoming paper [13]. 

4. On-line simulation. There are many ways we could make a multifrac- 
tal time-change of a CEBP. However, by defining the time-change via the 
crossing tree, we obtain a fast on-line algorithm to simulate the process. 
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As before, we will suppose that we are given subcrossing orientation distri- 
butions p^ and weight distributions F^ a , satisfying Assumptions 2.1, 3.1 
and 3.2. Let Y be the corresponding MEBP. Then we will simulate the se- 
quence {(7?, Y(7~®))}- That is, we will simulate Y at the spatial scale of 1. 
Given the multifractal nature of the process, the choice spatial scale is not a 
restriction, as the process can be scaled to any desired resolution. An imme- 
diate consequence of the definition of the crossing times is the following 
bound on Y: 

Y(t) e (y(7?) - l,y(7?) + 1) for t G (7M+i)- 

The basis of our simulation is a Markov process, which describes the line of 
descent of the current level zero crossing, from the spine down to level 0. For 
n> m and k > let nim, n, k) be such that Cl n is a subcrossine; of C n , , , 

and let fl» e {1, . . . , ^ n+ljfc) } be the position of C n k within C^ n+l fc) . 

Using this notation, if n : i is the tree- index of CP, then for < m < n — 1, 
=S^ k) . Let = (^^fc),^,^^^, 

fe A which is a description of the level n super-crossing of C®, and 

the family it belongs to. 

Let N(k) be the smallest n such that k(0, n + 1, k) = 1, and put 

y(k) = (y°(k),...,y N ^(k)). 

Lemma 4.1. 3^ is a Markov process. 

Proof. We first show how to update y(k) to obtain y(k + 1). Let M 
be the largest m < iV(A;) such that 

S K(0,n,k) = Z K(0, n +l,k) for n = 0, . . . , m. 

That is, for all m < M we have that C m , n , n is the last level m crossing in 
its family. 

If M = N(k), then N(k + 1) = N(k) + 1, and y gains the component 
yW+VQc + l). Let n = N(k + l). Then we have k(0,ti, k + 1) = 2, ^ = 2 
and «(0,n + l,Jfe + 1) = 1. The distribution of (Z? +1 , A^ +1 , B% +1 ) depends 
on y(k) only through a" = ^4" +1 (l), which is given by A\{Z^). 

If M < iV(fc), then for n = M + 1 we have «(0,n,fc + 1) = k(0, n, fe) + 
!. ^(o,n, fc+ i) = ^(o,n,fc) + 1 and + 1,* + 1) = «(0,n + l,k). Thus 

K(0,n+l,fe+l) «(0,n+l,Jfc)' /c(0,n+l,fc+l) re(0,n+l,fc)' a U fc(0,n+l,fc+l) 

T->n+l 

^K(0,n+l,fc) ■ 

For n > M + 1, we have y n {k + 1) = ^"(A;). 

For n = M, . . . , 0, we generate 3^ n (fe + 1) recursively. We have n{0, n,k + 
1) = k(0, n, k) + 1, SJJ( 0in>fc+1 ) = 1, and «(0, n + 1, fc + 1) = k(0, n + 1, fc) + 1. 
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b y that 



The distribution of (Z' 




depends on y(k) only through ajgi +ljfc+1) = ^ { oi + 2,k + i)( S ^o,M+i,k+i)) = 



K(0,M+2,fc) v° K (o,M+l,fe)+l > ■ 
That y is Markov follows from the conditional independence of the (ZP, A2, 
Ku) given the orientations a£. □ 

From y(k), we get the orientation ajj of C2, and the weights 



To calculate the crossing duration Vf, we also need i?™ +1 (l), for n = 0, . . . , 
N(k) and W£. Keeping track of the spine weights i?" +1 (l) is no problem. 
Calculating is less straightforward. We do have that the Wj? are condi- 
tionally independent given the a® , but we do not have an explicit formulation 
of the density of (Wj?|a£ = ±). 

The simplest way to approximate the Wj? is to generate a BRW (using 
Pa and F R ,) for a fixed number of generations, m say, and sum the node 
weights across the final generation. However, this is exactly the same as 
setting the Wj? to be constant, then scaling the resulting process by 2~ m , 
so we will just set equal to its mean v a °. 

Remark 4.1. Writing Y as X o M ~ 1 , where X is the CEBP correspond- 
ing to Y, we note that X and A4 are, in general, dependent. However, in the 
case where X is Brownian motion, we can construct A4 independently of X , 
simply by taking the orientations a9 as i.i.d. random variables, equal to + 
and — with equal probability. This is because for Brownian motion X(TP) 
is just a simple random walk. In fact, in this case, there need not be any 
relation at all between the crossing tree of X and that used to construct Ad . 

4.1. Pseudo-code. We give pseudo code for simulating {(7^°, Y(7^))}, 
with the crossing durations T> n , \ approximated by p n , \K(W n . ;) (i.e., putting 
W n: i = v l , where i = a n: \). 

Updating y(k) is handled by procedures Expand and Increment. Proce- 
dure Expand checks if M = N(k). If so, it then generates the component 
y M+l (k) and updates N(k). Assuming M < N(k), procedure Increment 
updates y n (k) to y n (k + 1) recursively, for n = M + 1, . . . , 0. The actions of 
Expand and Simulate are illustrated in Figure 3. 

Given sample position Y(T£), sample time and crossing state y(k), the 
procedure Simulate applies the procedures Expand and Increment, calcu- 
lates Y(Tk + i), T^ +1 and y(k-\-l), then increments k. Procedure Initialize 
generates an initial Y(l~x), 7j° and y(l) suitable for passing to Simulate. 



■M+l 



R 



,n+l 

K(0,n+l,fc) 



(S n Ki o, n ,k)) forn = 0,...,iV(/c). 



22 



G. DECROUEZ AND O. D. JONES 



Level 3 




/ \ Level 2 

A, \ K ° Level 1 

// / \ \\ / I \ \)0 ■ Level 
Me 

FlG. 3. Action of the procedures Increment and Expand. Suppose that at iteration k we 
have generated the tree given by solid black lines only, so that N(k) — 1, and y N( - k '(k) 
describes the family of node ii. When we reach node ii, we are at the end of level 
and 1 crossings. To generate the next level node, we first need to increase N(k) by 
1 and generate the family of node 13, which is the role of the procedure Expand. Next, 
procedure Increment goes down the tree and generates the families of nodes '14, and '15, 
hence generating i§. 



Recall that u = P(a™ +1 = +K = +) and v = P(a™ +1 = +|a? = -). Here 

N(k)+1 . . , A N{k)+l f y,N(k)+U 
a l is given by A 1 l [Z 1 ). 

Procedure Expand y{k) 

s N(k) _ z N{k)+l 

11 D K{0,N(k),k) ~ K(0,7V(fc)+l,fc) inen 

K(0,N(k) + 2,k) = 1 

Generate a^ k ^ +2 using u, v and a^ k ^ +1 

Generate (Z N{k)+2 A N{k)+2 E N ( k )+ 2 ) 

venerate ^ re ( 0j jv(fe)+2,fe) ' A K,(0,N(k)+2,k) ' a K(o,N(k)+2,k) > 

using the distributions p l A and 

conditioned on the first offspring having orientation 

where i = a^ k ^ +2 £ {+, — } 
q N(k)+i _ , 

J K(0,N(k)+l,k) 1 

Store < (fc)+2 (l) 
N(k) = N(k) + l 
End If 
End Procedure 

Procedure Increment y n {k) 

7^ Assume that C^ -1 is at the end of a level n crossing, 

# so S "(o!n-i,fc) = Z K(o,n, k y This is alwa y g the case for n = 
k(0, n,k + l)= k(0, n,k) + l 
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If S« = Z n ,y Then 
re(0,n,fcj K(0,n+l,fe) 

Increment X n+1 (k) 

on _ i 

J K(0,n,Hl) 

Generate (-£™(o, n +i,fc+i)> Ac(o,n+l,fc+i)' ^«(o,n+i,fc+i)) 
using the distributions and F^, 

Where i = <J 2 n + 2, fc+ l)( S "( + o'n + l,A: + l)) ^ {+, "} 

Else 

Xi(k + l)=yi(k) for g = n + l,...,JV(fc) 

en _ on , -i 

°K(0,n,fc+l) °K(0,n,fc) ~ 

End If 
End Procedure 

We apply procedure Increment to y°(k), and then it is recursively applied 
to all y n (k) such that C^ ? is at the end of a level g + 1 crossing for all 
< q < n. y n (k + 1) = y n {k) for all n larger than this. 

Procedure Simulate 
Expand y{k) 
Increment X°(k) 

Put i = ^(0,1^+1) ('S'jfc+i) 
If i = + Then 

Y(T k ° +1 ) = Y(T k °) + l 
Else 

End If 

7? +1 = T k ° + t,' nS +1) (<oW + D ( 5 
fe <- fc + 1 
End Procedure 

To initialize the algorithm, the procedure Initialize is used. Recall that 
(v + ,v~) T is the right /u(l)-eigenvector of M(l). 

Procedure Initialize 3^(1) 

k = l, N(1)=0, «(0,0,1) = 1, k(0,1,1) = 1 
Put a\ = i = + with probability a 
Generate {Z\,A\,R\) using the distributions 
p l A and FL, with i = a\ 
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Store R\(l) 
T? = v i 

If i = + Then Y(7f) = 1 Else Y(lf) = -1 End If 
End Procedure 

An implementation is available from the web page of Jones [16]. An ex- 
ample of the type of signal obtained with this algorithm is given in Figure 4, 
where we have represented an MEBP process with its corresponding CEBP. 
and F^ a are described in the caption. 

140 | , , , , 1 




-20 1 ' ' ' 1 1 1 1 ' 1 

0.5 1 1.5 2 2.5 3 3.5 4 

x10 4 



Fig. 4. Top figure: CEBP process where the offspring consist of a geometric (0.6) number 
of excursions, each up-down or down-up with equal probability, followed by either an up-up 
or down-down direct crossing [compare this with Brownian motion, for which there are 
a geometric (0.5) number of excursions]. Bottom figure: MEBP process obtained from a 
multifractal time change of the top CEBP process, with i.i.d. gamma distributed weights. 
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4.2. Efficiency. Consider the tree descending from crossing C 1 down 

to level 0. On average has fj, N ( k ) level subcrossings, so we must 

have N(k) = O(logfc). At each step, the number of operations required by 
procedure Expand is fixed [independent of N(k)], but we can go through 
Increment up to N(k) times, so the number of operations required by 
Simulation is of order N(k). Thus, to generate n steps, we use 0(n log n) 
operations, since ^J! =1 log /c = 0(n log n), and O(logn) storage. The algo- 
rithm is on-line, meaning that given the current state [of size O(logn)] we 
can generate the next immediately [using 0(logn) operations]. 

5. Randomizing the starting point. Crossing times are points where the 
behavior of the process can change, spatially and temporally, and the higher 
the level, the more dramatic this can be. For MEBP processes, is a crossing 
time for all levels, and because of this we cannot expect MEBP to have 
stationary increments. To avoid the problem of being special, we would 
like to start the process at a "random" time, as if the process had been 
running since time immemorial and we just happened across it. 

To make the idea of a "random" starting time more precise, let Y be an 
MEBP and {Y^} the nested sequence of processes used to construct Y, 
where Y^ is a single level n crossing from to ±2 n . Choose a time t 
uniformly in [0,T™] = [O,V n]0 ]. For any i 6 T n:0 , the probability that t is 
in C n -i is proportional to the crossing duration T> n . \ = p n .;VV n: i. That is, 
choosing t is equivalent to choosing n:j G dT n : so that the probability 
that n:j||;| =n:i is proportional to Pn-.iWn-.i- It turns out that we can do 
exactly this using a size-biased measure for a multitype branching random 
walk. 

Size-biased measures for branching processes were introduced by Lyons, 
Pemantle and Peres [27] and generalized to branching random walks by 
Lyons [26]. Kyprianou and Sani [21] then extended their construction to 
multitype branching random walks. Fix n, and for brevity write i for n : i. 
Let Q be the space of marked trees, where the mark associated with node 
i is (— logi?i| fc _ 1 (i[fc]), aj), writing k for |i|. Let T be the <7-field generated 
by all finite truncations of trees. The offspring orientation distributions 
and weight distributions F R , induce a measure £ on (Q, T). Let (} be the 
space of trees with a distinguished line of descent i £ <9Y, called a spine, 
and J- the cr-field generated by all finite truncations of trees with spines. 
Kyprianou and Sani define a size-biased measure fr on (Cl, J 7 ) such that 

(5.i) ^ (a( r,j) = ^ m . 

This is precisely what we want, and, remarkably, the measure can be con- 
structed using the original multitype branching walk, modified so that the 
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offspring generation down the spine is size-biased. That is, rather than con- 
struct y( n ) and then choose a spine, we can construct the process and the 
spine together. 

Let x 6 <9Y be the spine, and let p A and F^, be the offspring orientation 
and weight distributions for nodes on the spine. Then from [21], Section 2, 
we have that 



\a\ 

r- 



7=1 Js ^ r 



Note here that s and r are in Putting r = oo' a ' to get p l A (a), and then 
dividing out p A (a) to get F|j| a (r), gives us 



i«i , 

3=1 Js ~ r 

That these are well defined follows from Assumption 3.1. 

In the case where the offspring weights are i.i.d. with distribution F, we 
get 

f A (a) oc \a\p A (a), 



*fl|»«E / sF(ds)l[F(r(i)) 



The first of these is clearly a size-biased version of pV The second can be 
interpreted as conditioning on which offspring is on the spine, then size- 
biasing the weight for that offspring. 

For selecting the next node on the spine, we again have from [21], Sec- 
tion 2, that 

p a , r (j) := P(x[n + 1] = j\A x]n = a, R^\ n = r) oc v a ^r(j). 

Kyprianou and Sani also also show that under fr, the sequence {a x |J™ =1 
of orientations down the spine is Markovian, with transition probabilities 

v+ ° y 1 Mm ( v+ ° 



MULTIFRACTAL PROCESS WITH EMBEDDED BRANCHING PROCESS 27 



The stationary distribution is (u + v + , u v ), and so the reversed chain (mov- 
ing up the spine) has transition matrix 

^ {1 „°-)~W(" + --)' 

and the same stationary distribution as before. Note that it follows from 
assumptions 2.1 and 3.1 that u, v > 0. 

5.1. MEBP construction with random start. We now show how, given an 
MEBP Y : [0, oo) — > M generated by and F R ,, we can construct a shifted 

version, Y: (—00,00) — >M, with a "randomly" chosen starting point. Where 
unambiguous, we will use the same notation to describe Y as Y, and we 
will assume that assumptions 2.1 and 3.1 hold throughout. As before, we 
start by constructing a crossing of size 1 (level 0). Let x be the spine, which 
will be the line of descent corresponding to time 0. Accordingly, we will 
write Cq U = C x | n for the level — n spinal crossing. Note that previously, the 
first crossing at level — n was labeled 1, and started at time 0. For our new 
construction, time will occur somewhere in the interior of crossing Cj~ n , 
so crossing C-f n will still be the first full crossing to occur after time 0. 

The generation n (level —n) nodes in T n are totally ordered according to 
the rule i < j if and only if, for some m, i| m =j| m and i[m + 1] < j[m+ 1]. 
For i,j G T n let 

r |{k:i<k<j}|, i<j, 
d(i,j) = 0, i=j, 
[-|{k:i>k>j}|, i>j. 

We will write CZ n for Q. 

d(x,i) 1 

Set the orientation of Cq to be + with probability u + v + , and then gener- 
ate (Aq, Rq) using p A and F l R \ a , where i = olq. Choose j G {1, . . . , Z } using 
V A m , Rzi an d then put x|i = j. Subsequent generations are produced using 
and Fft< a for nodes off the spine, and and F R , for the spinal node. The 
spinal node in the next generation is chosen using p ar . Crossing durations 
are defined as before; that is, T>^ n = p^ n yV^ n , where W\ is the 7r-a.s. limit 
of ^jex„nT- Pi/ Pi- F° r ^ 7^ (nodes off the spine) the convergence of this 
sequence a.s. and in mean follows as before. For k = (nodes on the spine) 
a.s. convergence follows from (5.1) and the fact that p x \ n W x \ n £ (0, 00) £-a.s. 

Given crossing durations, we define crossing times as follows. Time 
corresponds to the spine x. For any m > 0, T§ m > is the first time the 
process starts a level — m crossing: 

V m = lim Yl 

ieT„,i| m =x| m ,i>x 
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rr—m 

>k+l 


q — m 
_ ' k 






= lim 

n— >oo 








ieT n ,i m =x| m ,i<x 


q — m 

'-k-1 


q — m 

— '-k 


- pZT^Zk for k > 1. 



We also put Y(0) = and 

y(V m )=iim_ £ «i r 



n— >oo 

ieT„,i| m =x| m ,i>x 

^(7^7) = y(T fc " m ) + a^2~ m for fc > 0, 

n— »oo * — ' 

ieT„,i| m =x| m ,i<x 

^(Tr^x) = ^(7r fc m ) + al^2" m for fc > 1. 
So for k>l, C fc " m is from T fc ~ m to T k ~^, while for A; < it is from 

to tt™. 

Let y (0) be the level crossing constructed above. We now show how to 
extend the construction from Y ^ to Y ( n+1 ) . Let n:x be the spine start- 
ing at level n. First choose ag +1 = i using the reversed Markov chain 5.2, 
then choose (Aq +1 , Rq +1 ) and (n + l):x[l] = j using p^, and p ajJ ., all 
conditioned on ckq, which is the orientation of (n + l):x[l]. Put the jth 
level n subcrossing of y( n+1 ), that is Cq , equal to Y^ n \ For the other level 
n subcrossings, we use the construction of Section 3.3, and scale the feth 
subcrossing by Rq +1 (k) / Rq +1 (j) . That is, we use the weights up the spine, 
from level to n, to rescale the process. Let Y be the limit of the fW. 

To see that Y(t) is defined for all t £ M we need two things. First we note 
that from the form of p a ,r, with probability 1 we cannot have n:x[l] equal 
to 1 eventually, or equal to Z n - eventually. That is, at all levels there will 
be crossings to the left and right of the spinal crossing. Second, we need to 
know that the scaling coming from the spine weights grows to infinity, that 
is, n£.Ro(A::x[l]) ^0 a.s. asn^oo. 

As noted above, the sequence of orientations up the spine is a Markov 
process. Because the weights are conditionally independent given the orien- 
tations, the sequence (^^Jlogi?Q(A;:x[l]),aQ +1 ) is Markov additive. Thus, 

Y2t\ log-R§(fc : x[l]) -» -oo a.s., equivalently njS Ro(k ■ x[l]) -> a.s., pro- 
vided the expected increments of the sum are negative. That is, provided 
the following assumption holds (this replaces Assumption 3.2). 



Assumption 5.1. Let i? ± be a random spinal weight, chosen according 
p A , F^ a and p a ,r- Then we assume that 

u + v + ElogR + + u'v'ElogR" < 0. 
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It remains an open problem to show that the process Y has stationary 
increments. 

5.2. On-line simulation. To simulate Y we need only modify procedures 
Expand and Initialie. Note that the spinal crossings are now counted as 
crossing at each level, so N(k) is the smallest n such that k(0, n + 1, k) = 0. 

Procedure Expand y(k) 

Whlle S K (0,N(k),k) = Z K (0,N(k)+l,k) D0 

K(0,N(k) +2,k) = 

Generate a^^ +2 using (u + t> + , u~t; _ ) and a^^ +1 
Generate Aq^ +2 , Rq^ +2 and Sq 1 
using the distributions p^, and p a ^ T 

conditioned on offspring Sq^ +1 having orientation a$^ +1 

where i = cXq^ +2 € {+, — } 

Store Jtf « + a (5 f 
N(k)=N(k) + l 
End While 
End Procedure 
Procedure Initialize 3^(0) 

k = 0, N(0) = 0, k(0, 0, 0) = 0, as(0, 1, 0) = 



Put ak = + with probability u" 1 



v 



Generate A$, Rq and Sq using the distributions 

Store i^(S§) 
7? = 0, F(7?) = 
End Procedure 
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